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Abstract 



Heavy Baryon Chiral Perturbation Theory (HBChPT) to leading order pro- 
vides a kernel to solve the Bethe-Salpeter equation for the P33 (A(1232)-channel) 
tt — N system, in the infinite nucleon mass limit. Crossed Born terms include, 
when iterated within the Bethe-Salpeter equation, both all one- and some two-pion 
intermediate states, hence preserving elastic unitarity below the two-pion produc- 
tion threshold. This suggests searching for a solution with the help of dispersion 
relations and suitable subtraction constants, when all in-elasticities are explicitly ne- 
glected. The solution allows for a successful description of the experimental phase 
shift from threshold up to y/s = 1500 MeV in terms of four subtraction constants. 
Next-to-leading order HBChPT calculations are also used to estimate the unknown 
subtraction constants which appear in the solution. Large discrepancies are encoun- 
tered which can be traced to the slow convergence rate of HBChPT. 
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1 Introduction 



The study of meson-baryon scattering is one of the most challenging topics of particle 
physics at low and intermediate energies. Firstly, because there exist accurate experimen- 
tal data and secondly because it provides an excellent scenario where Chiral Symmetry 
Breaking (ChSB) can be tested. Unlike the light pseudo-scalar meson-meson scattering, 
the role played by ChSB in the pseudo-scalar meson-baryon scattering case turns out to 
be more subtle and the construction of an Effective Field Theory (EFT) to describe the 
latter processes becomes more cumbersome. In the light meson-meson sector the expan- 
sion parameter is p 2 / (47r/) 2 , being p M the four- momentum of the Goldstone bosons [|lj and 
/ the corresponding meson weak decay constant. Since the baryon mass is already of the 
order of 47r/, such an expansion does not work in practice for the meson-baryon sector, 
though some efforts were undertaken 0. This drawback was overcome by formulating a 
new consistent derivative expansion for baryons in a chiral effective field theory || ||, || , 
called Heavy Baryon Chiral Perturbation Theory (HBChPT). HBChPT is a reformula- 
tion of Chiral Perturbation Theory (ChPT) with baryons where the baryon mass is shifted 
from the propagator to the interaction couplings of the EFT Q In this way higher orders 
are suppressed either by powers of p 2 / (4irf) 2 or by powers of p/M, where M is the baryon 
mass and p is either the light pseudo-scalar meson mass or the off-shellness of the baryon, 
respectively. Recently, a fully relativistic formulation consistent with the power count- 
ing proposed in HBChPT has been proposed ||, but there is so far a lack of practical 
applications. 

Standard HBChPT has been applied to the ttN system several times already in subse- 
quently higher orders || ^|, [10], [Tl|] and successful fits in the close to threshold region have 



been found. On the other hand HBChPT, as it is also the case for ChPT in the meson 
sector, is unable to describe resonances, which clearly indicate the appearance of non- 
perturbative chiral dynamics. Resonances, play a crucial role in meson-baryon systems. 
For instance, the low and intermediate energy irN scattering is mostly dominated by the 
A(1232) resonance, since it is close to threshold and its coupling to pions and nucleons 
turns out to be quite strong fl2"fl . Hence, the A degrees of freedom are often explicitly 
incorporated in any effective approach aiming at describing low and intermediate energy 
pion-nucleon dynamics to some degree of accuracy J7j-fl3|. 



Exact Unitarity plays a crucial role in the description of resonances, since at the cor- 
responding energy the scattering amplitude takes its maximum possible value. HBChPT 
only restores unitarity perturbatively and hence even the computation of phase-shifts 
becomes ambiguous beyond perturbation theory. Actually, in previous works the unita- 
rization program of HBChPT for the lowest S and P partial wave amplitudes obtained in 
has already been studied up to third order in the chiral expansion, by using the Inverse 
Amplitude Method (IAM) either in its most straightforward version ]TJ| or within an 



improved IAM ||15|f |. We note here that since the unitarity correction takes first place at 

3 This program was firstly suggested, consistently formulated and experimentally tested in the case of 
heavy quark (Q = b, c) physics [fj. 

4 The literature on unitarization methods is vast. References in the present context can be traced from 

Ref. ca. 
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Figure 1: Diagrammatic representation of the BSE equation. It is also sketched the used kinematics. 



NNLO in HBChPT one needs to go to that order in order to unitarize a la IAM. 

The aim of the present work is to show how a unitarization framework based on 
the Bethe-Salpeter Equation (BSE) JTB] and HBChPT to lowest order yields to a very 
simple expression for the scattering amplitude which preserves unitarity and provides 
a satisfactory description of tcN elastic scattering data. Our work can be regarded as 
a natural extension of the approach previously proposed in Refs. [f[7] and |IS| for tttt 
scattering to the physics of the irN system. 

We have concentrated our attention into the P33 {L-^nj) channel, and show that the 
existence of the A resonance can be understood in terms of the irN interactions con- 
strained to HBChPT. Hence there is no particular need to include the A-degrees of free- 
dom explicitly in order to provide a successful description of the experimentally observed 
phase-shifts. The extension to other partial waves as well as higher order contributions 



in HBChPT will be analyzed and presented elsewhere [19 . 



2 The Bethe-Salpeter Equation and the Elastic P33 
ttN Amplitude in HBChPT. 

The BSE for the elastic ttN scattering amplitude (T 1 ) in a given isospin channel reads 
t P (q,q') = v P (q,q')+iJ^/ P (k,q')D(k)S(P-k)v P (q,k) (1) 



with 



T P (q,q')^i = u a ,(P-q')t p (q,q') Ua (P-q) 

ViWW = u a ,(P-q')v P (q,q')u a (P-q) (2) 

where the kinematical variables are given in Fig. [IJ and V 1 , D and S are the two particle 
irreducible amputated Green function, pion and nucleon propagators respectively, t 1 , v 1 
and u are matrices and a spinor in the nucleon Dirac space and finally, a, a 1 are nucleon 
spin indices (helicity, covariant spin, etc.). The normalization of the amplitude T 1 is 
determined by its relation with the center or mass (CM) differential cross section, in the 
isospin channel /, and it is given by 



da 

— [q',a <- q,a 



(20 V 1 ' J 64vr 2 s 



Tp(q,q)<ra' (3) 



5 We use the normalization uu — 2M. 
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where s = P 2 . Rotational, parity and time reversal invariances ensures for the on shell 
particles 

2p(g, g'W' = "Svrv^ {^(s, 6)5 aa , + i B\s, 6) (n ■ a) a , a } (4) 

with 9 the CM angle between the initial and final pion three momentum and n a unit 
three- vector orthogonal to q and q'. Partial waves, f/ L (s), are related to A, B by 



A 1 = E 

L 



Pl(cos6) 

dP L (cos9) 
d6 



(5) 



The phase of the amplitude T 1 is such that the relation between the partial wave ampli- 
tudes and phase-shifts is the usual one, ff L (s) = e 1 s il( s j sin5j L (s)/\ q\cM- Hence, on-shell 
unitarity implies for s > (m + M) 2 

t r^j/M 1 I _^ I \^(s,M 2 ,m 2 ) 

MfiM]- 1 = -I q\cM = (6) 

with M = 938.27 MeV, m = 139.57 MeV, the nucleon and pion masses respectively, and 
X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 

Note that, to solve Eq. (TJ) both the off-shell potential and amplitude are required. 
Clearly, for the exact potential V 1 and propagators D and S, the BSE provides an exact 
solution of the scattering amplitude T 1 [|TIJ. Obviously an exact solution for T 1 is not 
accessible, since V 1 , D and S are not exactly known. We propose an expansion along the 
lines of HBChPT both for the exact potential (V 1 ) and the exact propagators. 

The lowest order of the HBChPT expansion in the P33 channel leads to || [1C[| 

i=— 1 i=— f i=— t 

Vp 2 (?) Q )a,a' ~ Vp 2 (q, q ) CT)CT / (crossed + Vp 2 (q, q (contact 
Vp~ 2 (q,q')a,a' {crossed = 2 ( ^ J U a >{Mv + p) ^ ^ ^ ? . U a (Mv + p) 

\f J v-(p-q')+ie 

i=- 1 / 

V P 2 (q,q')a,a' (contact = -j^ u a , (Mv + p') (v ■ q + v ■ q') u a (Mv + p) 

D{k) » l -— Sip) » l — (7) 

k l — m z + le v ■ p + ie 

where the velocity v is a time-like unit four-vector, in terms of which the nucleon momenta 
can be written as P — q (') = Mv + p (') with v ■ p (') << M . The EFT can be written in 
terms of nucleon fields N v (x) with a definite velocity t> M which are related to the original 
nucleon fields ty(x) by N v (x) = e lM ^ v x ^/(x). On the other hand, = ^ r y50' t j,uV 1 ' is defined 
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in terms of the velocity and the Dirac matrices^. Finally, / = 92.4 MeV and qa = 1-26 
are the pion weak decay constant and the vector axial coupling, respectively. Note that 
because of isospin conservation, the direct Born term does not contribute to the I = § 
potential. Likewise due to angular momentum conservation the contact term does not 
contribute either to the P33 ttN scattering and we will ignore it from now on. Thus at 
lowest order, the BSE (Eq. ([I])) for the P 33 channel reads: 

1 1 S-qS-k (g) 



k 2 — m? + ie v ■ P — M — v ■ k + ie v ■ P — M — v ■ (k + q) +ie 

We look for solutions to the above equation, suggested by a one-loop calculation, of the 
form 

7=2 

tp 2 (q,q) = A(v ■ q,v ■ q ) S ■ q S ■ q + B(v ■ q,v ■ q ) S ■ q S ■ q (10) 

where A, B are functions of two variables to be determined. Note that, as a simple one 
loop calculation shows, there appears a new dependence (S ■ q' S ■ q) not present in the 
lowest order potential V I= ?. That is similar to what happens in standard HBChPT || pTof] . 
The Eq. (Q) determines the functions A, B, which turn out to satisfy the following integral 
equations 

'g A \ 2 1 1 (9aV . f d 4 k k 2 -(vkf 



^' g ' W) = 2 7 v.P-M-v.(a + a>) + ie-3[T) 'I 



f J v ■ P - M - v ■ (q + q') + ie 3\f J J (2tt) 4 k 2 - m 2 + ie 
1 A{y ■ k,v ■ q') 



v-P-M-v-k + ie v ■ P - M - v ■ (k + q) + ie 
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B(v ■ q,v ■ q' 




d 4 k k 2 -(vkf 



(2tt) 4 k 2 - m 2 + ie v ■ P - M - v ■ k + ie 

A(v ■ k,v ■ q') - Bjv -k,v -q') 
X v • P - M - v ■ (k + q) + ie [ ' 

The above set of linearly coupled integral equations needs to be regularized, solved and 
renormalized. This program is highly non-trivial, and hence to draw solutions to these 
equations we will make use of the analytical structure of the involved functions. 



3 The following properties are satisfied Q| 



3 1 

S-v = S 2 =-- {Sfj,,S u } = —- {g^v - v,iV v ) [S fl ,Sv]=ie fl „ pa v p S a (8) 



with e i23 = +1 • 
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Figure 2: Diagrammatic representation of the BSE equation obtained by iterating only the crossed term 
potential, defined in Eq. (Q). Vertical lines represent the irN — ► irN and nN — * 7T7t7V cuts. 

3 Elastic Unitarity and On-Shell Dispersion Rela- 
tions. 

For on-shell scattering and at leading order in the 1/M expansion, implicit in HBChPT 
formalism, LAB and CM systems coincide. Then, choosing = (1,0) and taking into 



account that in the infinitely heavy nucleon limit v-q = v-q' = v-P — M = J m 2 + q CM = 
(s — M 2 + m 2 )/2y/s = d, the pion energy, we find 

/£M = -l( W 2 -mVH (13) 

where A(u) = A(u,u) for on-shell scattering. The function B enters in the P31 pion- 
nucleon scattering amplitude, and we will discuss it elsewhere |BJ. Elastic unitarity, 
Eq. (H), applied to the function A(u) requires for 10 > m 

lmA^(uj) = —(uj 2 - m 2 )f (14) 
v > 24tt V ; v ; 

From Eq. (|TT|), the on-shell function A(oS) satisfies, for u > m, 

d A k k 2 -(vk) 2 



A(u>) 




(2tt) 4 k 2 - m 2 + ie 



X 7—- V (15) 

c<j — u • /c + ie v ■ k - 



The function A (a;) above develops an imaginary part from two difference sources: the 
elastic channel and the inelastic channel corresponding to tcN — > 7T7riV process, as can 
be seen in Fig. (Q). To obtain, Eq. (|L3|) we have completely neglected all possible in- 
elasticities. We are going to proof that the solution of the Eq. flT5p automatically satisfies 
the elastic unitarity condition, given in Eq. fli~4]) , provided the pion energy, u, is such that 
we are below the two pion threshold or, equivalently, if the latter contribution to ImA(u) 
is neglected. This can be easily seen by standard operator methods by writing Eq. (|TT|) 
in an operator form 

A{u,uj') = U(u,u')+ f duU(u,u)G(u)A(u,u') (16) 
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with an obvious identification of the kernel U and G{v) given by 

^ * r dAk k 2 -{vkf i 

= -6 7 W V-r* + **-v-k + * * {V - V ' k) (17) 

with a — v ■ P — M. If the inelastic 7riV — > ttttN channel is not considered, the potential 
U{oj,uj') can be considered to be real, so that the corresponding ie can be ignored. Thus, 
one gets 

A(u,u')-A*(u',w) = I dvA*{u',v)2ilm.G(v)A(u,v) (18) 



By direct application of Cutkosky's rules, we get from Eq. fllTD 

2iImGH = -- / \k 2 -(v-k) 2 } {-2m) 2 5 + {k 2 -m 2 )5{a-vk)5{v-v-k) 

6 J (27r) 4 l j 

= -(m 2 - a 2 )5(a - i/)2iImJ (a) (19) 
6 

where # + (A; 2 — m 2 ) = 8(k°)5(k 2 — m 2 ) and the one-loop function 

f d 4 k 1 1 

Jo(a) " "V (2ic)*k*-m* + iea-vk + ie (20) 

has been defined. It is linearly divergent and hence two subtractions are needed to make 
it convergent 0, 

J (a) = K + K ia + J (a) (21) 



— \ Op 1 — TTIj Oi 

Jo(a) = 7, {arcosh i7r}; a > m (22) 

4n z m 

where Kq and K\ are related to the divergent constants Jq(0) and Jq{0). Finally, in the 
infinite nucleon mass limit, uj = uj' = a, and Eq.(^8|) reduces to the unitarity condition of 

An exact solution of the BSE (JTSJ) contains both elastic and two pion inelastic con- 
tributions, even though Eq. ([15|) constitutes a lowest order approximation. This feature 
makes a practical solution of that equation very difficult, mainly due to the non- locality 
in the momentum variable of the exchanged nucleon. On the other hand, we have seen 
that the solution of Eq. (|15|) exactly satisfies the unitarity condition of Eq. fll4]) when in- 
elasticities are neglected. Besides, the inelastic contributions start at (9(1// 4 ), but they 
are not complete because terms of the same order have not even been included in the 
iterated potential and propagator (vertex correction and self-energy insertion). Thus, we 
propose to solve the BSE (|15D in the limit in which the inelastic ttN — > ttttN channel 
does not contribute to the imaginary part of A(u), although this channel may contribute 
to the real part of A(u) as a subthreshold effect. To do that, we make use of the analyt- 
ical structure of A^iuo) in the complex u plane implied by the BSE to the lowest order 
considered. Once that the two pion channel contribution to lmA(u) is neglected, A~ l (uo) 
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only has a right hand cut, with branch points in u = m and infinity , and a discontinuity 
through the cut determined by Eq. fll4j) , 

A~ l {oj + ie) - A-\u - ie) = 2iImA(u; + ie) = — (u 2 - m 2 )i (23) 

On the other hand A~ l {uo) is not expected to have poles which would correspond to zerosQ 
of A(u) and thus of the cross section. Under these circumstances and as demanded by the 
asymtotic behaviour deduced from Eq. (|23|) , , a four-times subtracted dispersion relation 
for yl -1 (co>) is needed, which then readsQ for u > m: 

A-\u) = ^ + P( w ) + ( w 2 - m 2 )J ( w )/6 
z 9a 

PH = m 3 (co + Cl (--l) + c 2 (--l) 2 + C3(--l) 3 ) (24) 
V m m my 

The coefficients Cj might be expanded in powers of 1/f 2 . In the former equation we have 
explicitly separated the lowest Born term [—uf 2 /2g\) which only affects to the linear 
coefficient in u. Thus, we have 00,1,2,3 = 0(1) in the infinite nucleon mass limit. We 
would like to make a few remarks: 



The structure of the solution, Eq. fl24"|), of Eq. ( [T5D , agrees with the findings of 



Ref. 18 , where it was shown that in the case of ChPT in the meson-meson sector, 



to get the on-shell amplitude, the off-shellness of the BSE could be ignored, as long 



as a renormalized potential is iterated. With this philosophy in mind, Eq. (|i~5|) 
becomes 

A' 1 (to) = H(u) <r 



2g A 



^ ~ 6 J (2tt) 4 k 2 - m 2 + ie u - v ■ k + ie ( ' 

The divergent Integral H(u) is related to Jo(cu), H(u) = {uj 2 — m 2 ) J (u>) / 6 + Pi(u) , 
with P\{ui) a polynomial in the variable u of degree one. The above equation, is 
equivalent to Eq. (p4"|). Note, that in this approximation, the nucleon propagator 
which enters in the iterated potential, see Eq. (^), has been taken out of the inte- 
gral and hence it can never be put on the mass-shell. This is the reason why the 
approximate solution given in Eq. (p5|) does not contain the inelastic channel Nun 
contribution (see Fig. 0) and thus it can coincide with the solution of Eq. (0), based 
only on the elastic unitarity requirement. 

If one is going to neglect explicitly the inelastic channel, then only the nucleon 
propagator — v ■ k + ie) contributes to the elastic cut. This suggests to make 



7 Note that the kinematical zero at uj — m has been factorized out 

8 Note that, in this case, v4 _1 ([j) and (cu 2 — m 2 )Jo(w)/6 have the same analytical structure, and 
therefore they can only differ by a polynomial in w. 
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a Taylor expansion around uo in the integral of Eq. (|T5|) , it is to say 

k 2 -{vk) 2 A{vk,u) k 2 -u 2 A(u 



k 2 — m 2 + ie v ■ k — ie k 2 — m 2 + ie 

The leading order of this expansion yields again to Eq. 



+ 0(v -k—u) 



(26) 



4 Numerical Results 



To prove the usefulness of our main result, Eq. (|24f) , we perform a direct x 2 fit to the 
experimental phase shift in the P33 channel |20| from threshold up to ^/s =1500 MeV. 
We have assigned a 3% uncertainty to the phase shift as in [T(J plus a systematic error of 



one degree to the data (similar treatments are followed in [13| and [15 



The x 2 fit yields the following numerical values for the parameters: 



c = 0.045 ±0.021 



ci = 0.29 ±0.08 



c 2 = -0.17 ±0.09 



c 3 = 0.16 ±0.03, 



(27) 



with x 2 /d.o.f. = 0.2. The description of the data from threshold to the region well above 
the resonance is pretty good as can be clearly seen in Fig. ^ (upper solid line). 

On the other hand, the above coefficients can also be obtained from the 0(1/ f 4 ") 
HBChPT pieces obtained in Refs. || and [ID"|. Matching the BSE scattering amplitude 
to that obtained in either of the above references, and in the infinite nucleon mass limit, 
we find 



67rm 2 u; 2 



or 



m 2 )g 



x (mt 



(3,3) 



4:WJ (uj){t 



(1.1)12 



(2* 



Pol. of degree 3 in [u>— m) 



where the dimensionless amplitudes t^ 3 ' 3 ^ and t^ 1 ' 1 ) are defined in Ref. |15|]. They depend 
only on u/m, Qa and the Low Energy Constants (LEC's):6i + b 2 , &i6 — &15 and 619 in the 
notation of Ref. . The right hand side in Eq. (ffif ) is not a polynomial by itself, because 
it contains chiral logs stemming from the left hand cut, and it has to be Taylor expanded 
around u = m to comply with the polynomial structure of Eq. (|24|) . 



Using the values obtained in Ref. |TD| for the LEC's and translated into the notation 
of |Q as shown in Table I of Ref. |L5| (Set II entry) we find after Montecarlo propagation 
of errors, the lower solid and dash-dotted lines in Fig. |3| and estimates for the coefficients 

C 0,l,2,3 



,th 



0.001 ±0.003 



Jth 



0.038 ±0.006 



,th 



0.064 ±0.005 



„th 



0.036 ±0.002 



(29) 
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Figure 3: P33 phase shifts as a function of the total CM energy yfs. The upper solid line represents a 
X 2 -fit of the parameters of Eq. (|24| ) to data of Ref. |2(J (circles). Best fit parameters are given in Eq. (|27|). 
The lower lines stand for the results obtained with the parameters given in Eq. (^jj). Central values lead 
to the solid line, whereas the errors on Eq. (p9|) lead to the dash-dotted lines. 



Obviously, the structure of the Eq. (^1) which ensures exact elastic unitarity, allows 
for a satisfactory fit to data, as long as the contributions coming from two pion production 
are small and the left-hand-cut contribution can be modeled, in the scattering region, by 
a polynomial of the third degree in the variable uj — in. 

In the extreme static approximation, the HBChPT estimate of the polynomial P(u>) 
provides the bulk of the resonant shape, although it is far from describing accurately the 
data. Similar conclusions, within the IAM scheme, were reached in Ref. There also, 
it was shown that the 1/Mf 2 corrections (NLO in the chiral expansion) are extremely 
important. This might explain the big discrepancies between the set of parameters of 



Eqs. (|27|) and ( p9|) . Results in Eq. (p7|) are obtained from a best fit to data and effectively 
incorporate, buried in the fitted parameters, all higher order corrections, in particular the 
above mentioned 1/Mf 2 pieces. 



5 Conclusions 

In this paper we have dealt with the study of the A(1232) resonance as it arises in 
ttN scattering in the P33 channel. We have used HBChPT as a perturbative guide to 
implement in a sensible way the constraints imposed by chiral symmetry. The Bethe- 
Salpeter equation provides, in addition, a workable scheme where unitarity is exactly 
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restored up to the level imposed by the corresponding kernel and, moreover, allows for an 
identification of the infinite set of diagrams being summed up. In the case of 7rN scattering 
in the P33 channel we have used HBChPT to lowest order as input. The BSE is able, out 
of its divergence structure, to generate some of the higher order counter-terms required 
in standard HBChPT in the form of subtraction constants for the divergent integrals. 
Actually, our BSE amplitude admits a very good fit to the experimental phase shifts in 
terms of four subtraction constants. Nevertheless, their fitted values differ, at least within 
estimated errors, from those deduced from HBChPT. The discrepancy may be understood 
because the NLO (1/Mf 2 ) contributions are comparable to the LO (l// 2 ) ones, and hence 
the matching procedure of the BSE solution to standard HBChPT amplitudes induces 
large numerical corrections. This situation is not new and was already encountered in 
previous studies based on the IAM [|14], [15] which require a similar matching procedure of 
unitarized amplitudes to the standard ones of HBChPT. As we see, however, the amount 
of work required to describe the irN scattering data within the BSE, and in particular, 
the A(1232) resonance is drastically smaller than that needed for the IAM. In addition, 
as we pointed out in Ref. [Tj| for the meson-meson problem, one can search for systematic 
improvement and convergence in a BSE framework. A full study of the ttN system along 
these lines is left for future research mM. 
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